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We give a detailed calculation for the Hamiltonian of single galileon field theory, keeping track 
of all the surface terms. We calculate the energy of static, spherically symmetric configuration of 
the single galileon field at cubic order coupled to a point-source and show that the 2-branches of 
the solution possess energy of equal magnitude and opposite sign, the sign of which is determined 
by the coefficient of the kinetic term 02- Moreover the energy is regularized in the short distance 
(ultra-violet) regime by the dominant cubic term even though the source is divergent at the origin. 
Wo argue that the origin of the negativity is due to the ghost-like modes in the corresponding 
branch in the presence of the point source. This seems to be a non-linear manifestation of the ghost 
instability. 



I. INTRODUCTION 

Recently there has been a flurry of interest in the gahleon modification of gravity[l][2][3]. The theory 
is motivated by the DGP branc vi^orld model [4]. Galileon fields can be considered to be a generalization 
of the boundary effective field theory on the DGP brane at the so called Decoupling limit [19] where it 
corresponds to the brane-bending mode (for a review see [7]). The simplest version of galileon fields are 
postulated to have a novel symmetry structure, ic the action is invariant under tt — > tt + a^a;^ + b (a^, b arc 
constants). Demanding that the equation of motion only contains derivatives of order up to 2, gives rise to 
unique Lagrangian terms at each order in tt up to total derivatives and arbitrary coefficients. Remarkably 
the highest order, n, is determined by the number of dimensions, d, of the space-time, where n = d+1. This 
theory has been extended to include several independent galileon fields [13] [12] [5] [17] [2 7], supersymmetry 
[20], curved space generalisations[30][31] and covariant completion[28][29]. Although motivated by the DGP 
model, the theory is interesting and peculiar in its own right. There are novel field theoretic properties 
both at the classical and quantum mechanical level [24] [13]. In particular it is possible to choose suitable 
parameters to avoid ghost instabilities in the self-accelerating branch as opposed to the DGP model where 
there is no freedom to choose these parameters appropriately [1]. Furthermore violations of null-energy 
condition can be obtained without causing any instability [25] [26]. 

In this paper we derive the Hamiltonian for a single galileon field living in Minkowski background space- 
time with an arbitrary time-like boundary at spatial infinity. This has previously been done for multi-galileons 
without taking into account the boundary contribution [5]. Here we keep careful track of all the boundary 
terms and investigate the energy of the static spherically symmetric galileon field at cubic order sourced by 
a point-mass at the origin. We find that the energies for the non-trivial and normal (Minkowski) branch 
have equal magnitude but opposite signs depending on the sign of the coefficient of the quadratic term 0:2 
(see (2)). Setting a2 > gives positive (negative) energy for the normal (non-trivial) branch and vice versa, 
indicating ghost like behaviour in the branch with negative energy as we discuss later. This is a non-linear 
manifestation of the perturbativc ghost instability that has been explored extensively. 

Section-1 illustrates the framework used in computing the Hamiltonian. We use the normalization of 
energy with respect to a reference solution as was done in [10] [18]. In section-2 we present the ADM 3-1-1 
splitting for the bulk Lagrangian density. In scction-3 we do a subsequent decomposition of the boundary 
terms. We present the general expression for the Hamiltonian in section-4 and in the final section we use 
this result for a static spherically symmetric single galileon field and explore the implications of this result. 
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A. Infrared regularization of the Hamiltonian 

Our aim is to calculate the Hamiltonian for single galileon field theory living in Minkowski space-time 
with closed boundary(scc figl). The boundary is made up of constant-time hypcrsurfaces at far-past and 
far-future, S_oo) ^-i-oo and bounded by an arbitrary time-like hypersurface, B, at spatial infinity, with no 
inner boundaries. Usually it is fairly straightforward to calculate the Hamiltonian from the action of a field 
theory, where the Hamiltonian is the Legendre transformation of the Lagrangian, but it is slightly non- 
trivial when the action has boundary terms as in GR (Gibbons-Hawking- York boundary term). We follow 
a method that is conceptually similar to that followed by [10] in defining a physically meaningful notion of 
Hamiltonian for unbounded space-times. This is done by regularizing the action with respect to a reference 
field as explained below. 

The most general action for a single galileon field, Tr{x), in 4-D is given by[l][ll], 

^galileon ~ ^bulk ~t~ ^boundary (1) 

where, 

n— 5 



Sbulk = / (^) 

i" = {-a„7r„,7r[«^<3. ..<:!} (3) 

n— 5 p 

Sboundary = 2)77^ 7753 TT Trfj . . . TT^' | (4) 

Here tTo = dair and n±, tTq^ are orthogonal and tangential derivatives with respect to the boundary. We 
use the convention that antisymmetrization over the a indices do not involve the prefactor ^. Note that in 
the subsequent sections index, a, runs over 0..4 and i runs over 1..3. This is a consistent action for dirichlet 
boundary condition, ie when the fields and their tangential derivatives are held fixed at the boundary [11]. 
Sboundary is the analog of the Gibbons-Hawking- York boundary term in GR. Note that a2 > as we have 
defined in (3) yields a stable Minkowski branch free of ghost-like behaviour due to the positivity of the 
kinetic term, however this would make some other branches unstable. A concrete example of this is infact 
what we discuss in the final section. The action defined above is finite for compact geometries but diverges 
for non-compact space-times. To renormalize this action for non-compact space-times we choose a reference 
background ttq that asymptotes to the value of tt and also a solution of the theory. Then we demand the 
physical action to be given by, 

^physical — Sgalileon ['^] SgalUeon [^o] (5) 

consequently the physical Hamiltonian is, 

^physical ~ ^galileon [^] ^galileon [^o] (6) 

In order to derive the Hamiltonian for galileon field theory, one must do an ADM decomposition of the 
action. We postpone the final result until we have presented the decompositon of the bulk-space-time and 
the decomposition of boundary terms in terms of relevant derivatives and geometrical quantities. 

B. Bulk Decomposition 



We start with decomposing the action for single galileon field in terms of time and spatial derivatives. 
This is similar to the ADM formalism in GR except there is a prefered time direction since we are working 
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in Minkowski space-time. We consider the galileon field in Minkowski space-time bounded by a time-like 
boundary at spatial infinity B (see FIG. 1). We foliate the bulk space-time in constant-time space- like 
hypersurfaces E^. Thus the natural embedding is as follows, 

i:t:[x'']^[t,x'] (7) 

where t, define the standard cartesian coordinates giving the line element as, 

ds^ = -dt^ + Sijdx'dx^ (8) 

Ignoring the boundary terms (4), the most general Lagrangian for galileon fields in 4-dimensions can be 
expressed as follows, 

5 

Lgalileon = (9) 

n=2 

We consider a general term L" of order n in tt and seek to do a 3+1 split in terms of time and space. In 
the spirit of integrating by parts, we rewrite the Lagrangian as a piece that contains no 2"'* order time 
derivatives, L^uH,, and a total derivative term, -^^J^/t-o^er' (^^® Appendix [I] for details). Thus, 



^ — ^bulk "T ^left-over 



(10) 



where, 



here,"C2 



r(") in/~i -2 [is ii in] \U ia »r»] 1 



(")(»-!) 
2 



^7eft-over ~ 



{n-2){n+l) 



di. 



• 2 H-i 14 i 



+ {n- 2)d' 



7r,7ri3 7r'^7r>...7ri;; 



- (n - 2)8" 

+ {n - 2){n - 2>)d,^ 



7r„7ri3 7r^^7r.*...7r.^ 



(11) 
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Inserting the boundary term (4) back into the action and using Stoke's Theorem to convert bulk-integrals 
to boundary-integrals and including terms of all-orcler in 4D we recast the total action as follows, 



where 



^ total — ^bulk ^total — boundary 
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Sbulk -Y^htl «n{"C2 7T27r[>-....7r;:l -7r,,7r[-7rl^..7r;:l} 

n=2"' "'^t 



(12) 



(13) 



'^total-boundary — ^ ^ ^total-boundary 
n=3 



with, 



on 

^ total — boundary 



an j dt j 



(n-2)(n + l) 



77 TT^ TTj.TT 



(n — 2)7rr 



(14) 



(15) 



+ (n - 2)(n - 3)7r7ri,ri3 7r'*3 7r*^7rJ^^ . . . < 
+ a„ / I (n - 2)7ry Iw^^^wlt • • • - ^i^^^'Kl ■ ■ ■ | 

JdM L J 

Note that Sf^f^,_i^^^^^^^y = 0. Here Try = y^aTr, tt^ — r"9a7r denote the derivatives along the normal 
vectors y^r" (see FIG. 1) respectively. tTq denotes the covariant derivative with respect to the boundary, 
B (see next section). Also, — J dtN Jg^ where A'^ = (1 + 6'^)~^^'^,6 = Uav"' is the lapse function. We 
have completed the decomposition of the bulk-terms in the action. In the next section we decompose the 
boundary term Stotai-boundary with respect to the relevant derivatives to be defined below. 



C. Boundary decomposition 

We seek to decompose the boundary terms in terms of derivatives with respect to the closed 2-surface 
St = B f\Y,t and derivatives along r"',U°'. We work in full-space time coordinates and begin by presenting 
the definitions of various derivatives and projection operators. 

Tab = gab + UaUb '■= Projection operator forS^ (16) 
Hab = 9ab — VaVb ■= Projection operator fori? (17) 
Qab = Hab + riaUb = 'Jab - TaU ■= Projection operator for 5* 

We use Da, Da, Da to denote covariant-derivatives with respect to T,t,B,St- For brevity this convention is 
used on the indices in long expressions. 7r„, Try, TTrj^TT are derivatives along the corresponding vector fields 
defined as Tr„ := D„Tr := n°_Da7r etc. Also, Tr„a := DaDnir, TTrd ■= DaDri^. 7r„2 := D^^Tr := DnDnir, Tr,.2 := 
D^n := DrDrTT. The action of a covariant derivative Da on a hypersurface (with an associated projection 
tensor hab) on a given tensor lying on the surface is given by [22], 



DoT^^:::^; = hX[ ...hXhi\... h^brt::i; (is) 

Boundary terms contain derivatives Da, Da, DaDj,, DaDb which can be decomposed as follows (see Appendix 
[11]). 
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D^-K = DaTT + raD^TT (19) 
Dair = DaTT - UaDnTT 

DaDbTT = DaDbTT - Klf,DnTT - 2n(„£)6)£>„7r + 2rn^aKl)^D''TT + naUbDl-K 

Here if^j, K'^f^ = OK^f^ are extrinsic curvatures of the 2-surface St with respect to the hypersurfaces Sj, B 

respectively. T^ab) = '^{Tab + Tba) denotes symmetrization of the indices. We can now express the boimdary 
terms given in (16) by substituting the decomposition given above. Thus the boundary terms are (here we 
omit the result for 5*'' order for brevity, see Appendix [II]), 

Sltal-boundary = / _^ { " 3(1 + 0^)7r\r ' ' 3^(1 + e^f^TTTT^, - 6(1 + e^f'^TT^ + [D-Kf-Kr^ 

(20) 



^total-boundary — OH j ^ | - 2Tr^TT^°'Kl 

-2(l + ^2^-5(^7r)2 



(21) 



e^TT^ + 20{1 + 9^)hTrr + (1 + e^)7Tr^ 
0(1 + 6|2)7T7f + 26*2(1 + 0'^)h7rr + 0(1 + 6'2)7T7r^2 + (1 + ^2)1 TT^Tf 

+ 20(1 + 02)^^7r^ + (1 + 02)i7r^7rH - K^n {^^^^^ + (1 + 0^)4 + 2^(1 + 0^)hTrr) 

0{1 + 02)5,r27ra^" + (1 + 0^)TTr^naK + 2^(1 + 0')5i7r,7Ta< + (1 + 2^2)^^^^.^^ ^ 02^2^.^a 



lob 



e'^n^ + (1 + 0'')TT'f + 20(1 + 0")5^7r, 



q2\_2 



- 2(£)2^) 



^(l + r)57i-^ + (1 + 6* )7rr'' + (1 + S0^)n^TTr + 36'(1 + r)57T7r^ 



+ 2(1 + 02) 5(^^+(l + 02-)i^^-)^ 



0(1 + 02)'7r^+302(l + 02),rV+30^(l + 0^)57r/^+(0* - l)7r/ 



- bir'^TTr^K' - 67i-^(D^7r)7r^ - bn' K^t,TT"'TT" + 47T^7r™7r" - 47r/7r™7r° + 47r^2^ifc7r"7r'' 
+ 27rr3(-£,2^^ _^ 27r/i4:^ + 27r^7ra7r["7r?l + 27r^2^a7rl"i4:"^ + 27r^7r^2 (L»7r)2 - 27r7r^(L)7r)2 



D. Derivation of the Hamiltonian 

Having recast the galileon action interms of ADM decomposition we can now write down the Hamiltonian 
directly. The Hamiltoninan density of the galileon theory described by the Lagrangian L{'jt, dir, ddw) is given 
by the legendre transform. 



H = pn — L 

where the canonincal momenta p is given by. 



(22) 



6 



P 

I 111 i~ the Hamiltonian of single galileon field theory is, 

5 

Hgalileon 



dL 

dn 



5 

n=2 ''^t 



(23) 



(24) 



total — boundary ^ total — boundary total — boundary 

Where the last 3 boundary-terms are given by (20), (21) and (56) (Appendix [II]). 



E. Energy of static galileon fields coupled to a point-source 

Let us use our Hamiltonian to compute the energy of a single galileon field at cubic order in a static 
configuration with SO{3) symmetry, coupled to a point mass, m, at the origin. We take St to be a 2- 
sphere with fixed radius R. Here the theory contains two vacua: a normal branch (X^^) and a non-trivial 
branch {X^-^) (sec (28)). The stability of these branches depends on the sign of a2, where a2 > leads 
to a stable normal branch but an unstable non-trivial branch and vice versa. Here we demonstrate that 
this perturbative instability is consistent with our non-linear calculation using the full Hamiltonian where 
it manifests as negative energy for non-trival (normal) branch when a2 is positive (negative). The natural 
coordinates to work with are the spherical coordinates {r,9,(j>). The Hamiltonian function for this set-up 
becomes, 



H = A'K I drr'' { a27r'^ -|- 2a: 



(25) 



Here Mp is a dimension-full coupling constant with mass dimension, usually this is of order planck mass for 
gravitational theories. Also, (') = ^ and p = m5^^\r). Note that <S'^ota;-6oundary vanishes for this set up, 
since for static 5*0(3) symmetric galileon field, tt = DaTr = 0, and time invariance of the 3-boimdary, B, 
implies 9 := n.r = 0. The equation of motion is given by [l],(The it appearing in the expressions from (26) 
to (33) is the number tt not to be confused with the field.) 



a2X + SaaX^ 



(26) 



where X = —. The normal and non-trivial branch solutions of (26) are given implicitly by. 



r 



-a2 



3ma3 



6a3 



(27) 



-a2 



3ma3 



6a3 



The corresponding reference solutions which we choose to be the normal and non-trivial vacuum solutions 
are given by setting m = in (27). 



X 
X 







ref 
+ 

ref _ 0^2 
3Q!3 



(28) 
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It is convenient to rewrite the integrand in (25) using the equation of motion to ehminate the tt dependence. 
Thus, 



r m f 

H = -4tt / drr* {aX^ + 4a3X^} + dr {rX} 
Jo Mp Jq 

The energy for positive and negative branches is now given by, 



(29) 



E±=H[X±]-H[X''^f]\m=o 



Substituting (27), (28) above we get, 



135a§ ISMpQfa jo 



2MpTral 



27a^ 7o 



drr'^ 1 



3q:3T77, _3 



1/2 



(30) 



(31) 



After some change of variables the integrals can be recognized as a linear combination of hypergeometric 
functions given by. 



135q!^ 



+ 



3\/67ra3 



-1/2, 1/6, 7/6, ^ 



(32) 



— sign{a2) 



2od \ 27rm „7 



63 y 3Mpa| 



-1/2,7/6,13/6,-2^^"^^^" 



Sasm 



Here the hypergeometric functions are real and positive and defined for the range 0:3 > —' 3^ 
However for real values of as is forced to be positive. We now take the limit — )• 00 and the energy 

becomes. 



^00 



(33) 



We get a finite expression for energy with equal magnitiude and opposite sign. The infra-red divergence 
is regularized by substracting the vacuum energy contribution. As a non-trivial check for our calculation we 
take the limit m — >■ in (32) and obtain. 



lim E± = Q 

m— S-O 



(34) 



as expected. 



F. Discussion 



We conclude with a few remarks on our analysis of the energy of galileon field theory. Having presented the 
expression for Hamiltonian in ADM formalism carefully keeping track of all the boundary terms, we calculated 
the energy of static spherically symmetric configuration. In particular, the results of our calculation shows. 



• The two branches of the cubic theory coupled to a point source have energies of equal magnitude and 

opposite sign. 

• The expression for energy flips sign when the sign of 0.2 is changed. 

• Even though we couple galileon field to a divergent source at the origin, energy is still finite where 
non- linear cubic contribution dominates the divergent quadratic term and regularizes it. 

We argue that the negative energy of the non-trival(normal) branch when > 0(< 0) with a coupling 
to a point mass indicates a ghost like instability. Our calculations have been entirely classical and as was 
argued in [16] the appearance of negative energy can be traced back to the wrong sign in the propagator, at 
quantum level. If one evades negative probabilities by shifting the poles in the denominator of the propagator 
it leads to negative energy. Scattering processes involving ghost like particles and ordinary matter particles 
can generate ghost particles with unbounded negative energy and matter particles with unbounded positive 
energy. We believe the sign flip of the energy when changing the sign of a2 further reinforces this argument, 
for it is the correct sign of a2 in ordinary field theories that ensures the positivity of the kinetic term in 
the Lagrangian. It is interesting to note that a similar calculation was done for Gauss-Bonnet gravity in 
[18] and the authors found that the energies for the 2-branches match both in magnitude and sign. Further 
more it was shown that one of the vacua of Gauss-Bonnet gravity was unstable despite the fact that ghost 
like modes were not excited by the spherically symmetric black- hole [6]. In contrast here we find that point 
source which can be taken to be a spherically symmetric source in the limiting case, does seem to excite 
ghost-like modes giving negative energy. We would like to pursue this line of enquiry in future, it would be 
interesting to do this calculation for covariant galileon model and Multi-galileon theories. It is well known 
that the bulk part of the Hamiltonian vanishes identically for diffeomorphism invariant field theories [9] and 
it is not clear how this would play out for galileon models. In order to imderstand the origin of this negative 
energy it might be interesting to study the instanton transition amplitudes via bubble nucleation between 
the different vaccua in galileon models using the method pioneered by S.Coleman [8]. 



Acknowledgments 



I am deeply grateful to my supervisor Antonio Padilla without whose valuable and helpful suggestions 
this work would not have been possible. I would also like to thank him for reading through the manuscript 
carefully. 



G. Appendix I - Bulk Decomposition in Detail 

Consider the general n*'* order term appearing in the tt — Lagrangian. The highest order possible is (n+l) 
in space-time of n-dimensions. 

i^"^ = -«n{'ra.7r[«^C...<:l} (35) 

First we make note of the following general identities which we would make use of repeatedly. Note that 
Einstein-summation is assumed for repeated indices. 

aia2---an fi2---»n ' ' ' ' ' ili2---t...in ' ' ' ' ' ili2---in-lt ' ili2---in V^^/ 



t II 12 * 



zi \i i' 

rnr ^TT} TT,- 



(37) 
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Using (36) L" can be cast in the following form. 



L(") = -a„ {wtnKll . . . TT^l + {n- 2)1^,1,^^-4 ■ ■ • + ^i.^^'^^Z ■ ■ ■ } 



(39) 



Integrating by parts with respect to the upper time index in the second term before doing the anti- 
commutation operation we get, 



using (37) on the first term and (38) on the last term we get, 



(40) 



-Q!n|( 



(n - l)ntn'nl^ . . . - (n - l)(n - 2)^,^^^4y^^ . . . n^^ + w.^n^^-n^ • • • < 



+ {n- 2)9* 



(--2) 5.3 



(n-3)! 

again using (37) for the last term we get. 



= -a J {n - l)7rt7r*7ri^ . . . <1 - (n - l)(n - 2)7r*7r^3^[* _ _ _ ^ + _ 3)^* 



(»-2) 
(n-3)! 



t P4 » 



(n-2)(n-3) 
(r^-3)! 



a„<i C27r 7r]3 ... iTiJ - tt^^tt^ ^tt^ . . . tt^ 



^] (n-2)(n+l) 



• 2 fiQ 24 i,i 



(n - 2)5* 



(n-2)(n-3)9i; 



In the final step we have recast the following term as 



It lA I 



2 ^ ''[.^31 



TT, . . . TT. 

I4 l„ 



(41) 



and integrated by parts with respect is inside the commutator. Now we convert the term involving 9* [..] 
into a total derivative in full space-time by adding and substracting a corresponding term involving a total 
derivative with respect to the spatial slices S^. Thus we get, 



■ " __ 



(n-l)(n + l) 



da 



- (n - 2)8^" 
+ (n-2)(n-3)a„3 



'lju'laa 'I "04 • • • ''an 



+ {n- 2)8" 



2 

Ja3„a4 _a„ 



TT TT' -^TT^* •••^a„ 



(42) 



as promised. 
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H. Appendix II 

1. Decomposing the extrinsic curvature of B 

Extrinsic curvature of the time-like surface B is given by, 

K^b = K [d^^,] (43) 
we wish to decompose this interms of the following basis of one forms, 

Ev = Vadx" , Ea = Qabdx^ , En = Hadx" (44) 



we get, 



Thus, 



as expected. 



K^a = y^t^K^d = (45) 
K^v = V^VK^, = 

Kin = qln'Ki, = qln^Ht [daV,] = qtn^ [daV,] 

= qln'^Hi [dM] = qln" [dM] = V^n" [S^Kr"]] = (y.r)n'^ [d^ra] 



Kab = ^dS - 2 Hal if^g) + riaUbK^n (46) 



2. Decomposing the derivatives DaDtn, DaDbir 



First we derive the following results to be used later, 

Kb = QaDcTb = qllH [ddTe] = qH [ddVe] = qt [ddHre]] = qi [ddU] (47) 

Kl = ql [D,n,] = q^H^H^ [d,n,] = q^H^, [O^n,] = q'.q^ [dan,] = q^n, [OdM = {n.r)qi [d^r,] = {n.r)Kl, 

(48) 

D^DbTT = DM den] = qUbUql = qUt [dcddn] - qUt [ddVe] DrTT = qUl [dcdM - Kl,Drn (49) 
qSr'^ [dM = qtrb [dd{r''dc-K)] - qtrb [a^r^] d^n = [OaDrTr] - rbKl^D^n (50) 

raTbr'^r'^dcddiT = raVbr^ddr'^ddiT) - r aV braider '^)ddn = rarbD'^n - rarb{r''Dcr'^)dd'K = ranD'^n (51) 
we have used the result VaD'^Vb = in the last equality. 

qinbu'-'dcddTT = mqiddin^dcir) - nbq^{ddn'"}{dcTT) = nbDaDnir - n^gf [iJ^ + V''Ve]{ddn''){dcTT) (52) 
= HbDaDnTT - nbR'^^D^TT - nbqtYeddrf Dy-K = TibDaDnTT - HbRj^D'^'K + nbK?nDvn 

naUbTfn'^dcdd'n: = UaUbn'^dcin^ddTT) - nanb{rfdcn'^)ddn = UaUbDl-K - nanb[H^ + V'^Ve]{rfd^n'')dd-K (53) 
= naUbDn-K - nanb{rfDcn''')ddn - nanbU^Vedcn" DyiT = naiibDln + naUbRnnDyn 
we used naD^'Ub = Oin the last equality. 



11 



D.DtTT = [^^dcTT] = ^tlhdbM = iHd^da-K (54) 

= DaDbTT + + 2r(„|^|b)£),7r - 2r(„|i^|\),D% + r^n^l^ 

where(47),(49),(50), (51) was used. 

DaDt,n = D^[H-M = HiHlda[md,-K] = H^H^dcdan + HiHt{daHl){d,-K) (55) 
= qUbdcdd-K - 2qf^^n\h)n'dcddn + naUbrfn'^d.dd'K - H^H^dd[VeV'']dcW 
= ital^cddT^ - '^qfa\'^\b)rfdcdd'n + naUbrfn'^dcdd'n - K^b^V^ 
= DaDb-n- - BKlf^DniT - 2n(„|£)|6)£>„7r + 2n(„|/rp(,)^£)% + n„n6£»^7r 

where (45), (46), (47), (48), (49), (52), (53) was used. 

3. Boundary term at 5*^ order 



^total-boundary — ^^5 J '^^ j ^ | ~ ^'^^ 



(56) 



+ 3(1 + e'rHen + (1 + e')i7rr) - {nr^,^ + ^„7r[„|)7ri'^7r| - {0^7^^, + 4)KfKf 

+ 2i9\l - ^l)^,i,^K^^'Kf^ + {enl - 7r.^)7rjj<fl - 2(^7r^„a - ^W)Kfi,'^ 

- 2(^7r„ + nr)T^fKl''\,-K^ + {e\l - ■k^t.^-k^'^KI^KI'^^ + 27ra7r[«7rS7r„g + 27ra7r[»7r^]7r^g 

+ 29\a7r^''Kl'^Kl^TT'TT'^ + 2TiaTr^^Kf^Kln^n'^ ~ (On^ + 7rr)na7r^\^^Kl^^ - 2(7r„2 + ^,,2)^["7r^^Va 
+ 2(07r„7r„2 - TT^TT^O^aTrl^iff 1 - 2eTTlTTaTT^:Kl''^ + 247ran^^ K^"^ + 2{0^ttI - 4)t, aK^'' t,' 
+ 27r„7ra7r[f TT?' + 27r,7ra7r[«7rjl - 2(^7r„ + 7r,)7rf i^^^lTr^a 



+ 3 



+ TT.TTaTrl'^^^^Trf + ^^^aTrl"^,^^/^,^^! + 27r,7r^g7ra7r['^^^l - 2TTrTTa7r^^ K^"^ KI^w'tt'' + 2TTr7Tan^^Kl%\f_ 
+ 27r^ifi''7rg7ra7r[%^l + Tr^TTaTrl'^Trli^Jl + 27r^2 7r^7ra7r[«7r?l + 24nr2Tran^^ kI''^ - 2n^^naTT^^ kI''^ 

+ 24TTaKl^"Kl''K^ - 27r^7ra7r["7rgl + 2n^TTfKl"K^7ra 



+ 6 



If one needs to restrict the above expression to the basis -E„ = Uadx°',Er = radx°',Ea = qabdx^, the 
following expressions can be used to convert the relevant terms (we omit this step for brevity). 
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■Kn = {l + e'^)^TT + 6-Kr 

T^na = (1 + ^')'^a + 0T^ra ' K?„ [OtT + (1 + e^f-TT., 

7r„2 = (1 + e^)jT + 261(1 + e'^y^TTr + e'^'Kr2 - \{1 + 0'^)^ TTr + OtT 



(57) 
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